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I. INTRODUCTION AND SUMMARY 

Accelerated expansion of the Universe [1] is one of the 
most important discoveries in cosmology. Usually, it is 
explained by the existence of tiny cosmological constant 
or some scalar fields in the framework of the Einstein’s 
general relativity. However, it is worth exploring an al¬ 
ternative possibility that general relativity breaks down 
on cosmological scales and the accelerated expansion of 
the Universe results from modification of gravitational 
laws themselves. Many attempts have been made in this 
direction, but, it has been recognized that it is extremely 
difficult to construct a consistent theory. One of the prob¬ 
lems often encountered by the attempts to modify gravity 
in the infrared is appearance of ghost modes in the spec¬ 
trum of the theory. These are modes with the wrong sign 
in front of the kinetic term. They have negative kinetic 
energy and lead to vacuum instability with respect to 
creation of ghost particles together with ordinary matter 
particles, the interaction between ordinary matter and 
ghost being mediated (at least) by gravity. The rate 
of the development of such instability diverges due to 
infinite phase volume unless one introduces a Lorentz- 
violating cutoff, see e.g. [2]. 

An interesting model incorporating modification of 
gravitational laws at large distances was proposed by 
Dvali, Gabadadze and Porrati (DGP) in [3]. The model 
describes a brane with four-dimensional worldvolume, 
embedded into flat five-dimensional bulk. Ordinary mat¬ 
ter is supposed to be localized on the brane, while grav¬ 
ity can propagate in the bulk. A crucial ingredient of 
the model is the induced Einstein-Hilbert action on the 
brane. The action of the model is given by 

-I- J d'^Xy/^C 

mat •) 

where R is the five-dimensional Ricci curvature, 
is the intrinsic curvature of the brane computed using 
the induced metric 7^1,, and Lmat is the Lagrangian of 
the matter on the brane. Throughout the paper we use 


(—, -b, -b,...) convention for the signature of the metric. 
The induced gravity term is responsible for the recovery 
of four-dimensional Einstein gravity at moderate scales^, 
while at distances larger than 


Tc = 



( 2 ) 


gravity is five-dimensional. Cosmology in the DGP 
model is governed by the following modification of the 
Friedman equation [6], 
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where H is the Hubble parameter and Pmat is the mat¬ 
ter density on the brane. Two possible choices of sign in 
(3) give two branches of the cosmological evolution. The 
upper sign corresponds to the Universe whose expansion, 
in the absence of the cosmological constant on the brane, 
decelerates at late times, the Hubble parameter tending 
to zero as the matter on the brane dissolves. We call 
this branch of solutions the Friedman-Robertson-Walker 
(FRW) branch. On the other hand, the choice of the 
lower sign in (3) makes possible de Sitter expansion of the 
Universe with the Hubble parameter H = Ijrc even in 
the absence of matter. Thus, the latter branch contains 
the self-accelerating solution, where the accelerated ex¬ 
pansion of the Universe is realized without introduction 
of the cosmological constant on the brane. Below, we will 
refer to this branch of solutions as the self-accelerating 
branch. 

However, the self-accelerating branch of solutions turns 
out to be plagued by the ghost instability^. This was 
first demonstrated in [5, 7] using the boundary effective 
action formalism. In [8] this result was confirmed by ex¬ 
plicit calculation of the spectrum of linear perturbations 
in the five-dimensional framework. The analysis was per¬ 
formed in the case when the role of matter on the brane 


^ Let us mention that the mechanism of this recovery of four¬ 
dimensional gravity is rather non-trivial [4, 5]. 

^ This is not the case for the FRW branch. 
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is played by non-zero brane tension cr 7 ^ 0. At ct < 0 the 
ghost mode was identified with the scalar field describ¬ 
ing the brane bending, while at cr > 0 it was found that 
the ghost degree of freedom coincides with the helicity -0 
component of the graviton, which turns out to be mas¬ 
sive. By continuity, one expects the ghost to be present 
also in the case <7 = 0 which corresponds precisely to the 
self-accelerating cosmological evolution. 

The purpose of this paper is to show the presence of 
ghost in the self-accelerating Universe with cr = 0 by an 
explicit calculation. The analysis is subtler in this case 
than for cr 7 ^ 0. The reason is that at cr = 0 the masses of 
the graviton and the brane bending mode coincide and 
the two modes mix. As a consequence, it is impossible 
to diagonalize the quadratic Lagrangian for these modes 
and single out the Lagrangian for the ghost mode. In¬ 
stead, to demonstrate the existence of ghost, we use the 
Hamiltonian approach [9, 10]. We construct the Hamil¬ 
tonian for the helicity -0 excitations and find that it is 
unbounded from below. For modes of high momentum, 
k ^ H, the mixing terms in the Hamiltonian can be ne¬ 
glected and it decomposes into a sum of Hamiltonians for 
a positive energy mode and ghost. 

The paper is organized as follows. We start by consid¬ 
ering in Sec. H a simple model illustrating the subtleties 
encountered by the analysis of the spectrum in the self- 
accelerating Universe. In Sec. HI we find the pertur¬ 
bations about the self-accelerating solution of the DGP 
model and construct the four-dimensional effective action 
for the discrete graviton mode and the brane-bending 
mode. In Sec. IV we construct the Hamiltonian for the 
helicity -0 excitations and show that it contains a ghost 
mode. Technical details of the calculation of the Hamil¬ 
tonian are presented in Appendix. 


II. ABELIAN EXAMPLE 

To set the stage for the analysis of the quadratic theory 
for perturbations about the self-accelerating solution, let 
us consider a simple toy model whose spectrum exhibits 
analogous properties. We consider an Abelian gauge the¬ 
ory coupled to two Stiickelberg fields^. The Lagrangian 
of the model has the form, 

^ , (4) 

where /i, /| > 0. When the interaction between the vec¬ 
tor field and the scalar bosons is switched off by setting 
e = 0 , the fields cpi and (p 2 are an ordinary scalar field 
and ghost, respectively. Below we assume e 7 ^ 0. The 
system possesses gauge symmetry 

(p2 ‘P 2 + a , (5) 


which can be used to set ip 2 = 0- In the case ff 7 ^ /| 
one can diagonalize the Lagrangian (4) by introducing 
the vector field 


Bfi — A^ -\- 


if!-me ■ 


( 6 ) 


The result of the diagonalization has the form 
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fliy 


if? 



nf?-m 


. (7) 


It describes a massless scalar field and a massive vector 
with the mass = {ff — /|)e^. If /f — /| > 0 the 
scalar is ghost, while the vector field has positive mass 
squared. On the other hand, when ff — fi < 0 the sign in 
front of the scalar kinetic term is correct, while the mass 
squared of the vector becomes negative. This implies 
that the longitudinal component of the vector is ghost in 
this case. 

At the point /^ — /| = 0 in the parameter space the 
vector field becomes massless. Simultaneously, the trans¬ 
formation ( 6 ) becomes singular suggesting that the La¬ 
grangian cannot be diagonalized. In terms of fields A^, 
ipi the Lagrangian reads 

B - f?eAf,d^ipi . ( 8 ) 

The last term describes mixing between the vector and 
the scalar fields. Note that, if this last term were absent, 
the Lagrangian ( 8 ) would possess an Abelian gauge sym¬ 
metry which would make the helicity -0 component of the 
vector unphysical, and the theory would be ghost-free. 
However, the Lagrangian ( 8 ) as it stands, is not gauge 
invariant and the helicity -0 component of the vector field 
becomes a physical ghost through its mixing with the 
scalar field. 

To demonstrate this we consider the Hamiltonian of 
the theory ( 8 ). One concentrates on the helicity-0 sector 
as the transverse degrees of freedom are unaffected by 
the presence of the scalar field, and hence, their energy 
is positive definite. One performs the Fourier decompo¬ 
sition 


Aq 

Ai 


dfk 




(k,t) , 


(27r)3/2 

/ a(k t) 

J (2^)3/2 \k\ ^ ’ 


<Pk 

(27r)3/2' 


„ikx 




(9) 


Inserting these expressions into ( 8 ) and integrating out 
the variable n one obtains the Lagrangian for the fields 

Of, 


Lc\ — 






\m+m\k\acj>* 


( 10 ) 


^ We thank R. Rattazzi for suggesting this model to us. 
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where dot denotes derivative with respect to time. In 
writing (10) we made use of the relations q;(— k, t) = 
a* (k, t ), etc. Introducing the canonical momenta 


/le • _ /fe . .2 (. ihef 

““ \k\^' |fc| P 


we obtain the Hamiltonian, 

Ho = J + 




l ^|2 f2 


2 /?e 2 
— fickacj)* 


After the canonical transformation 

a = fiea, ^ = -fiea + fikcj) , 

'^a — fl^ i'^a ^</>) 5 — flklTfi^ , 

it takes the following form 

ay /’j3, + (/ie)^|~ |2 /y ^ 2 ~ 

Ho = / d-^ki ---Itt^I - [fie) 


-fc 2 + (/ie )2 
H-^- TTa 


( 11 ) 


( 12 ) 


(13) 


Clearly, this Hamiltonian is unbounded from below, neg¬ 
ative energies being associated with the field a. For the 
modes of high spatial momentum, k^ ^ (/ic)^, one can 
neglect the mixing terms in (13); then, it is clear that the 
mode a is ghost. 

At generic values of the momentum k, the Hamilto¬ 
nian (13) represents one of the normal forms of quadratic 
Hamiltonians. It cannot be diagonalized by a canonical 
transformation (see, e.g. [11]). This fact is a manifesta¬ 
tion of the resonance between the modes a and 4>. The 
solution of the equations of motion following from (13) 
has the form 










g±i|fc|t 


It contains a linearly growing part due to the resonance. 
Again, this part can be neglected for high frequency 
modes and at short time scales ^ ~ ^ ifieV ' 

solution in this regime becomes the sum of two purely 
oscillatory modes. 


III. PERTURBATIONS IN 
SELF-ACCELERATING UNIVERSE 

We now proceed to apply the Hamiltonian analysis, 
analogous to that of the previous section, to the DGP 
model. In this section we study the linear perturba¬ 
tions about the self-accelerating cosmological solution, 
and construct quadratic four-dimensional action for the 
localized modes. 


The five dimensional background metric corresponding 
to the self-accelerated branch has the form [ 6 ], 

= dy^ + N"^ [y)"fiii,{x)dx^dx'^ , (14) 

where N[y) = 1 -|- idjyj, and 7 ^,^ is the four-dimensional 
de Sitter (dS) metric with the Hubble parameter H. The 
brane is located at y = 0 , and Z 2 symmetry across the 
brane is imposed. The case of vanishing brane tension, 
which is of primary interest to us, corresponds to 



Tc 


However, we will not use this relation for somewhile in 
order to be able to compare the cases of vanishing and 
non-zero brane tensions. 

Let us now consider perturbations of the metric. We 
impose the gauge 


dgoo = Sgofi = 0, (15) 

and write g^^ = + hfii,[x,y). For perturbations 

obeying equations of motion it is possible to impose in 
addition to (15) the transverse-traceless (TT) gauge = 
0 , = 0 , where denotes covariant derivative with 

respect to the metric 7 ^ 1 /, and indices are raised using the 
metric . In this gauge the Einstein’s equations reduce 
to 

= 0 , (16) 

where □ = VaV'*'. We will call the coordinate frame 
where the metric is TT the bulk frame. In this frame 
the brane is displaced from the origin. We parametrize 
the perturbation of its position by y = ip{x). With the 
account for the brane bending one obtains the following 
junction condition at the brane in the TT gauge [ 8 ], 

— H [1 + Hrc)hf^i, 

= [2Hr,-l)[V^V^ + H^j^,)ip. 

Taking the trace of this equation results in the equation 
of motion for the field (p, 

Dip + iH^ip = 0 . (18) 

Let us summarize the results about solutions of Eqs. (16), 
(17) in the case Hvc 1 [ 8 ]. First, there are solu¬ 
tions leaving the brane at rest, tp = 0. They describe a 
tower of modes of the form h^i/(a;,y) = x^iw\x)Fm[y). 
The fields satisfy the four-dimensional equation 

= [mf + 2H'^)x^fw^ for massive spin -2 fields in 
dS space-time. The functions Fm[y) obey the following 
equations, 


„ vnf — 2H^ 


(19) 

F + F — 0 

~ ^ ’ 

y > 0 , 

JT, - 2HFm = -rn^TcFm , 

o' 

II 

( 20 ) 
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where prime denotes differentiation with respect to y. 
There is a continuum spectrum of modes with masses 
rn? > (9/4)i7^, and a normalizable mode 

FmAy) = (21) 


accelerating Universe we perform a field redefinition (cf. 
Eq. (6)) 

1 — 9 Hr 


with the mass 


ml = 


3Hrc - 1 


( 22 ) 


Second, there is a discrete mode with non-zero brane 
bending. The corresponding perturbation of the metric 
has the form, 


1 — 2Hr 

h^v{x,y) = ^ (V^V^ -|- H . (23) 

This is a solution with m? = 2H^. In [8] the effective 
quadratic action for the two discrete modes was con¬ 
structed for the case Hxc ^ 1- As the masses of the 
modes are different in this case, 7 ^ 2H^, the La- 
grangian decomposes into the Lagrangians for the mas¬ 
sive spin-2 field Xfj,u and the scalar field ip. It was found 
that at Hxc < 1 the scalar field is a ghost, while the 
tensor field is well behaved. [It is worth mentioning that 
the region Hvc < 1 corresponds to the unphysical case 
of negative brane tension. Nevertheless, considering this 
regime does make sense if the value of Hxc is close enough 
to unity.] On the other hand, at Hrc > 1 the ghost mode 
coincides with the helicity-0 component of the field 
while the scalar field has the correct sign in front of its 
kinetic term. At the point Htc = 1 the masses of the two 
modes coincide, and the solution (23) becomes singular. 
This signals that at the point Htc = 1 the Lagrangian 
for the discrete modes cannot be diagonalized. One con¬ 
cludes that the situation is analogous to the situation in 
the Abelian model of the previous section. 

Let us concentrate on the sector of the localized modes. 
Keeping these modes only, the expression for metric per¬ 
turbations at Hrc ^ 1 reads 


h^,.ix,y) =X^i^‘^\x)FmAy) 


1 _ OHr 

i?(r37^(v^v. + iLVMx). 


To obtain a non-singular Lagrangian for the localized 
modes in the limit Hrc t corresponding to the self- 


In terms of the fields (p the KK decomposition of 

the metric perturbation reads 


hfj,u{x.,y) — A.piiy{x)FYYi^{jj^ 

1 — 2Hr 

+ ^ + H 7^i/)(/3(x)(l — Fm^iy)) ■ 

Now, we can take the limit 1, keeping both fields 

A-^v{x) and ip{x) finite. Using Eq. (21) we obtain, 

= A^^{x) + log (1 -f Hy) . 

(24) 

We note in passing that this representation of h^i, in 
terms of localized modes can be obtained directly in the 
case Htc = t starting from the Einstein equations in the 
bulk and junction conditions on the brane. 

From now on we set Hvc = 1- The field ip still obeys 
Eq. (18). To obtain the equations of motion for the field 
we plug the expression (24) into Eqs. (16), (17). This 
yields 

= H (V^V, + ip . (25) 

Note, that the TT condition for the metric perturbation 
h^i, implies that A^^, must be also TT. 

We proceed to construct the effective quadratic action 
for the fields A^,^ and ip. The full quadratic action has 
the form: 

+ ^ j d'^x./^Sg‘“'SGi^i , (26) 

where are the variations of the five¬ 

dimensional and four-dimensional Einstein’s tensors, re¬ 
spectively. We impose the gauge (15) and write the ac¬ 
tion (26) in the Gaussian normal (GN) coordinate frame 
where the brane is at rest. 


J 


5(2) 



d^xV^(^ - + \hh" + 

+ J - Hh^ - ■ 


( 27 ) 
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Here h^ii, is the perturbation of the metric in the GN 
frame, h = h^, and 

X^,{h) = <5GW + 

= - i (□ V - V^Vah“ - V.Vo/i“ + V^V.h) 

- ]p^lv - Oh) + H'^ • 

In deriving the expression (27) we included the contribu¬ 
tions from the hrst integral in (26) which are proportional 
to 5{y) into the integral over the brane. Let us note that 
the action, restricted to the gauge (15) does not give (55) 
and (5/i) Einstein equations. So, a priori, the latter equa¬ 
tions must be imposed as constraints. However, we will 


see that the equations of motion coming from the effec¬ 
tive action imply that the metric perturbations are TT. 
This, in turn, implies that the (55) and (5/x) Einstein 
equations are automatically satisHed. 

The metrics in GN and bulk frames are related by the 
gauge transformation 

2N 

h^u = - 2HN-ffj,,yip . (28) 

To obtain the effective action for A^i,, ip one takes h^y in 
the form (24) and plugs the expression (28) in (27). Note 
that in spite of the logarithmic growth of the metric per¬ 
turbation (24) into the fifth dimension, the integrals over 
y in (27) are Hnite due to the presence of 1/N‘^{y) warp 
factors. After straightforward calculation one obtains 





H{A^^^Xf,Xy(p - AUip - SH^Ap) - , 


(29) 


where we introduced the notation A = A^.. Let us stress 
that it would be incorrect to impose the TT condition on 
the field in the action (29). In particular, assuming 
A/j^y to be TT would lead to the absence of the term 
describing mixing between the fields and p in (29), 
and hence, to incorrect Held equations. On the other 
hand, let us show that the equations of motion following 
from the action (29) entail both the TT conditions and 
the field equations (18), (25). By varying the action (29) 
we obtain 

- 2X^,(A) - 2H^A^y + 2H^^^yA 

- HX^XyP + H"f^yDp+ 3H^"f^yP = 0 , (30) 

Q/7 

- X^XyAf^^ -k DA -H 3H^A -(□ -k ‘lH^)p = 0 . 

( 31 ) 

Taking the covariant divergence and the trace of the first 
equation one obtains 

V^AO - V,A = 0 , (32) 

2Xf,XyA>^''-2nA + 3H{a + AH^)p = 0 , (33) 


where in deriving (32) we made use of the identity 
Xf^Xjf(A) = 0. Equations (32), (33) imply Eq. (18). 
Now, Eqs. (31), (32), (18) yield the TT conditions, 

V^A') = 0 , A = 0 . (34) 

With the use of Eq. (18) and the TT conditions equation 
(30) is reduced to Eq. (25). 

Before proceeding further, let us make the follow¬ 
ing comment. The first line in (29) coincides with the 
quadratic Lagrangian for the Fierz-Pauli theory of mas¬ 
sive graviton with the mass m? = 2H^ in dS background. 
Thus, it is invariant [10] under gauge transformations 

Afiy I ^ Afj^y + {X^Xy+H^^^y)coix) . (35) 

The mixing of the graviton field A^y with the scalar p 
in (29) breaks this symmetry explicitly. Let us restore 
the symmetry by introducing the Stiickelberg field tp. 
Namely, let us consider the following action, 
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S,_ff = 


5e// = ^ / + H^A^ 


+ 


■^ / r 1 2 

^([(V^V. + HV)'/^-— 


^(lU 


{D + 4H‘^)ip- ^A 


2H 




2H .i2 
3 

2H i2 


(36) 


(□ + 4i72)V,-^A]')| 


Note that the terms with four derivatives in (36) cancel 
out after integration by parts. The action (36) is invari¬ 
ant under the gauge transformation (35) supplemented 
by 


V2'-^'/5+—w, -Ip ^ ■i/’+—u; . (37) 

It is straightforward to check that in the gauge ip = 0 the 
action (36) reduces to (29). The form of the Lagrangian 
(36) is similar to the Lagrangian (4) of the Abelian model 
of Sec. II; it corresponds to equal coefficients in front 
of the Lagrangians for the Stiickelberg fields, /f = /|. 
In fact, one can check, that the quadratic action for the 
perturbations about the self-accelerating Universe can be 
written as an action for a massive graviton with the mass 
= 2H^ and two Stiickelberg fields in dS space-time 
at any value of the parameter Hvc- We will not pursue 
this issue further in this paper. 


IV. HAMILTONIAN FOR HELICITY-0 
EXCITATIONS 


The quantities in Eqs. (40) correspond to linear per¬ 
turbations of the spatial metric, gij = pSij + Aij, the 
shift functions, Ni = goi = A^i, and the lapse function, 
N = (—= 1 -I- n. The next step is to perform 
the spatial Fourier decomposition. We are interested in 
the helicity-0 sector, as we expect the ghost mode to be 
present there. (As to the helicity ±2 and ±1 sectors, they 
are not affected by the presence of the scalar if and do 
not contain instabilities [10].) So, we write 


= - 


d^k ,i,^f 
(27r)3/2 \ 



hkj\ 
1 ^ ) 


ai{k,t) 


kikj /, . 

<22 (k, t) 


(27r)3/2 \k\ 




(41a) 

(41b) 


< (41c) 

= J . (41d) 


Let us come back to the four-dimensional theory with 
quadratic action (29). Our aim now is to show that this 
theory has ghost in dS space-time. We adopt the ap¬ 
proach of Refs. [9, 10] and construct the Hamiltonian for 
the helicity-0 excitations. 

To simplify the formulas we set in this section the four¬ 
dimensional Planck mass equal to one, = 2^1 = 1. 
We use the following explicit form of dS metric: 

ds^ = —dt'^ + f^{t)6ijdx^dx^ , (38) 

where 


Note that we use the same notation n for the Fourier 
transform of the lapse function; this will not lead to con¬ 
fusion. The / factors on the l.h.s of Fqs. (41) are chosen 
in such a way that the time-dependence of the resulting 
Hamiltonian is simplified, see below. 

All but two degrees of freedom introduced in (41) are 
non-dynamical and are eliminated by making use of the 
constraints. The details of this procedure are presented 
in Appendix. Once this is done the Hamiltonian is ex¬ 
pressed in terms of two dynamical variables, the brane 
bending mode (p and the helicity-0 excitation of the gravi¬ 
ton V = ai — a 2 . The Hamiltonian has the form, 


fit) = e^‘ , (39) 

and i runs from 1 to 3. One decomposes various tensors 
into spatial and time components, and introduces 

Aij, Ni = Aoi, n = —Aoo/2 . (40) 

These notations are reminiscent of the standard notations 
in the Hamiltonian approach to gravity. Indeed, the Hrst 
term in (29) represents the quadratic action for gravita¬ 
tional perturbations in dS background with Sg^i, = 


Ho 



k* P 3\ 
3^4 + ^2 + 2 j 


2k^ 


l^v 


* 

4> 







(42) 


where the canonical momenta tTi,, are conjugate to the 
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variables f, (p, respectively, and 

■ (43) 

For all fields in the expression (42) one has f (—k, t) = 
u*(k, t), etc. Note that the Hamiltonian (42) explicitly 
depends on time via k. To simplify the Hamiltonian we 
perform a canonical transformation to variables v, TVy, 
etc. The transformation is conveniently parametrized by 
the generating function <P(v,(/>,7ry,n^;t) such that (see, 

e.g. [12]), 


9$ 


9$ 


etc. 


dnp ’ " dv* 

The Hamiltonian transforms in the following way, 

9$ 

Ho 1 -^ Ho + , 


(44) 


(45) 


where the partial time derivative acts on the explicit time 
dependence of the generating function. With the choice 


$ = 



6i72' 




+ 


2 


^vV* 

V6 


it^v* 

VQ 



(46) 


one obtains 


Ho = (rk 




- H’^TTyTTl 




dr 

2 


\3l 

2 


(47) 


From this expression, it immediately follows that the 
mode (j) is ghost. 

Two comments are in order. First, the Hamiltonian 
(47) is explicitly time-dependent via fc^(t) (see Eq. (43)), 
and hence is not conserved. Instead, the conserved en¬ 
ergy in dS space-time [9, 10] is 


E= = H+ / d^x T° Hxi , 


(48) 


where is the energy-momentum tensor of the system, 
and = {—l,Hxi) is the Killing vector of dS space-time. 
However, when one considers the system at distances 
much shorter than 1 /H (which is the case relevant to the 
ultraviolet stability), the second term in (48) can be ne¬ 
glected, and the energy coincides with the Hamiltonian. 
Second, the expression (47) is not diagonal, reflecting the 
resonance between the two modes. The situation again 
is similar to the case of the Abelian model considered in 
Sec.H. However, for the modes with wavelengths much 
smaller than the horizon size, ^ H^, one can neglect 
the 0{H^) terms in (47), and the two modes decouple. 
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APPENDIX A: ELIMINATION OF 
NON-DYNAMICAL VARIABLES 


Substitution of the Fourier decomposition (41) into the 
action (29) yields the Lagrangian for the helicity 0 sector: 


Lo = J - -jdi - 2Hai - 2Hn\^ - (di - 2Hai - 2Hn){a2 - 2i7a2 - 2\k\iy - 2Hn)* 

1 • ' ^ 

+ {(f>- H(j)){ai - 2Hai - 2Hn)* + ^{(t> - H(j)){a 2 - 2Ha2 - 2\k\iy - 2Hn)* - k'^pal + ylai]^ 


-b -b 2H‘^aial + H‘^\iy\‘^ + n* 


(aH^ + 2p) ai -b 2H^a2 - k^(, 


9 . 9^2 9ff2 


(Al) 



where k is defined by Eq. (43). In deriving (Al) we omitted total derivative terms and used the relations (ai(k))* = 
ai(—k), etc. The variable v can be integrated out explicitly. From the resulting Lagrangian one determines the 
canonical momenta: 


TTl = 


7r2 = 


'^(j) - 


dL 

da\ 

dL 

80*2 

dL 


= - 1 + 


IP 


3H 


3H 


di -—ai — 2Hn ] — { 82 -;^q ;2 — 2Hn + 1 + 




, ■ 3i? 

= — \ ai -~ 2 Hn 


- 


d(j)* 

and the Hamiltonian: 


3i7 

ai -Ofi — 2L/n 

2 


5H 

1 


3H 


7: a2 - - 2 Hn - - 


2 H^ 


5H 


5H 

-( 

2 


(A2a) 

(A2b) 

(A2c) 


'W frfij 1 I ^ V |2 1| |2 7 , 1 , 1 3H 3H 5H 

Ho = d k<- —IttiI + \j 2 ^ 1 P j “ 3''^'^' ~ ~ —7^20:2 + —’ 

+ k'^aicj)* - f ^ + hA |aiP - 2 H^aia *2 - + n* [ 2 Htti + 2 Htt 2 - + 2 k^)ai - 2 H^a 2 + 


Let us study the set of constraints of the system. From 
the Hamiltonian one obtains the primary constraint 


Xl = TTl + 772 - 


_l_ 2 H"^ 

— - ai - Ha 2 + —(j) = 0 . (A3) 


Taking the time derivative of Eq. (A3), using equations 
of motion and the relation 

we find the secondary constraint, 

TTl / 1 2^2 \ 2 4 2 

H H - 3“' - 3“^ = ° ■ 

(A5) 

Finally, vanishing of the time derivative of Eq. (A5) im¬ 
plies the constraint 


2n -I- 2ai -|- 0:2 = 0 . (A6) 

Note, that Eq. (A6) is equivalent to the tracelessness con¬ 
dition = 0. We use Eq. (A6) to eliminate the variable 
n from the Hamiltonian. The constraints (A3), (A5) are 
second class and have canonical Poisson bracket: 


{Xi,X2} = l- (A7) 

Let us reduce the number of degrees of freedom by per¬ 
forming a canonical transformation to variables tti, di, 
etc., such that 


^2 = Xi > 0:2 = X2 • (AS) 


To find a suitable transformation we introduce the gen¬ 
erating function <I> depending on the old coordinates 
ai,a 2 ,(t> and the new momenta ifi, 772 , such that. 


ai = 


d^ 


d<^ 


TTl = 


The requirements (AS) imply the following equations: 

5$ 54) 


7r2 = 

5$ 

dn^ 


da^ da 2 


fc2 -p 2 ff^ P 

-ai — Ha2 + 


H 


2 H 


1 5$ / 1 2 P 


H 


H dal 
2 54> 

3 

As a solution we choose 


54) 

^ ^ 


0.1 


- 0 L 2 ■ 


4> = d-^k 


2 P 

3IP 


771772 



9iJ3 


F2 


3H 


T^ 2 p(f, 


-\- 7T\0-y — 1^x02 ~\~ '^(pY 

o o 

P + 2 jj H, 3A:2 2 

4-—l«il + Haia2-—ai(j) 


The transformation generated by this function has the 
form. 


+ 0!i 0 : 2 , 

(AlOa) 

2 p f 1 2 p\ 2 


2 1 
+ gCH + 2«2 , 

(AlOb) 

4'=^i,+4.. 

(AlOc) 

2 p -i- Y p 

TTl =771 + 3^72 -f ^ ai -f Ha 2 ^^(j) , 

(AlOd) 

772 = - 77l -1- ^^2 + Hai , 

(AlOe) 


etc. 


(A9) 
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3k^ 


(AlOf) 


Using these expressions one can check explicitly that the 
conditions (A8) are satisfied. Making use of the con¬ 
straints we set 7r2 = 02 = 0, and we are left with two 
dynamical variables v = ai = ai — a 2 , 4> = 4> and 


their conjugate momenta. In the derivation of the re¬ 
sulting Hamiltonian one should account for the explicit 
time dependence of the generating function according to 
Eq. (45). A straightforward calculation yields the Hamil¬ 
tonian (42) of the main text. 
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